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The form of energy termed heat that typically derives from lattice vibrations, i.e. the phonons, 
is usually considered as waste energy and, moreover, deleterious to information processing. 
However, with this colloquium, we attempt to rebut this common view: By use of tailored 
models we demonstrate that phonons can be manipulated like electrons and photons can, thus 
enabling controlled heat transport. Moreover, we explain that phonons can be put to beneficial 
use to carry and process information. In a first part we present ways to control heat transport 
and how to process information for physical systems which are driven by a temperature bias. 
Particularly, we put forward the toolkit of familiar electronic analogs for exercising phononics; 
i.e. phononic devices which act as thermal diodes, thermal transistors, thermal logic gates and 
thermal memories, etc.. These concepts are then put to work to transport, control and rectify 
heat in physical realistic nanosystems by devising practical designs of hybrid nanostructures 
that permit the operation of functional phononic devices and, as well, report first experimental 
realizations. Next, we discuss yet richer possibilities to manipulate heat flow by use of time 
varying thermal bath temperatures or various other external fields. These give rise to a plenty 
of intriguing phononic nonequilibrium phenomena as for example the directed shuttling of heat, 
a geometrical phase induced heat pumping, or the phonon Hall effect, that all may find its way 
into operation with electronic analogs. 



PACS numbers: 44.10.+i, 63.22.-m, 66.70.-f, 05.70.Ln 
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I. INTRODUCTION 

When it comes to the task of transferring energy, na- 
ture has at its disposal tools such as electromagnetic ra- 
diation, conduction by electricity and also heat. The 
latter two tools play a dominant role from a techno- 
logical viewpoint. The conduction of heat and electric 
conduction are, however, two fundamental energy trans- 
port mechanisms of comparable importance, although 
never been treated equally in science. Modern infor- 
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mation processing rests on micro electronics, which af- 
ter the invention of the elect ronic solid-state transistor 
( Bardeen and Brattainl fl948h . and other related devices, 
sparked off an unparalleled technological development, 
the state of the art being electric integrated circuitry. 
Without any doubt, this technology markedly changed 
many aspects of our daily life. Unfortunately, a simi- 
lar technology which builds on electronic analogs via the 
constructive use of heat flow has not yet been realized by 
mankind, although several attempts have been repeat- 
edly undertaken. In everyday life, however, signals en- 
coded by heat prevail over those by electricity. There- 
fore, the potential of using heat control may result in 
an even more abundant and unforscen wealth of applica- 
tions. A legitimate question then is: Does phononics, i.e., 
the counterpart technology of electronics present only a 
dream? 

Admittedly, it indeed is substantially more difficult 
to control a priori the flow of heat in a solid than it 
is to control the flow of electrons. The source of this 
imbalance is that, unlike electrons, the carriers of heat 
- the phonons - are quasi-particles in the form of en- 
ergy bundles that possess neither a bare mass nor a bare 
charge. Although isolated phonons by itself do not in- 
fluence each other, interactions involving phonons be- 
come of importance in the presence of condensed phases. 
Some examples that come to mind are phonon polari- 
tons, i.e. the interaction of optical phonons with infrared 
photons, the generic phonon-electron interactions occur- 
ring in metals and semi-metal structures, phonon-spin 
interactions, or the phonon-phonon interaction in pres- 
ence of nonlinearity. Therefore, heat flow features aspects 
which in many ways are distinct from charge and mat- 
ter flow. Nonetheless, there occur in condensed phases 
many interesting cross-interplays as it is encoded with 
the reciprocal relations of the Onsager form fo r mass , 
heat and charge flow, of which th ermoelectricity ( Calient 
119601: iDubi and Di Ventrat l201ll ) or ther mophoresis, i.e . 
the Soret vs. its reciprocal Dufour effect ( Calient [l960h . 
are typical exemplars. Therefore, capitalizing on the 
rich physical diversities involving phonon transport as 
obtained with recent successes in nano-technology may 
open the door to turn phononics from a dream into a 
reality. 

With this Colloquium we shall focus on two fundamen- 
tal issues of phononics; - i.e. the manipulation of heat 
energy flow on the nanoscale and the objective of pro- 
cessing information by utilizing phonons. More precisely, 
we shall investigate the possibilities to devise elementary 
building blocks for doing phononics; namely we study 
the conceptual realization and its possible operation of 
a thermal diode which rectifies heat current, a thermal 
transistor that is capable to switch and amplify heat flow 
and last but not least a thermal memory device. 

The objective of controlling heat flow on the nanoscale 
necessarily rests on the microscopic laws that govern 
heat conduction, stability aspects or thermometric is- 
sues. For these latter themes the literature provides 



the readers with several comprehensive reviews and fea- 
tures already. For heat flow and/or related thermo- 
electric phenomena on the micro-/nano- and molecular 
scales we refer the readers to the excellent treatises b y 

Gafperin et all (120071). 



Casati (2005) 



Lepri et al. (12003). .. _ _ 

Dharl (l2008h iLi et all (l2 005rj).IZhang and Lil (l2010D.|Pop| 
(|2010l ). and recently also bv lDubi and Di Ventral (|201lh 
We further demarcate our presentation from the sub- 
jects of refrigeration o n mes oscopic scales and thermom- 
etry (|Giazotto et all , l2006f ) and, as well, as our title 
suggests, also do not address per se in greater detail 
the issue of conventional way of manipulating heat flow 
upon changing thermal conductivity by means of various 
phono n scat t ering mechan i sms in nano- and heterostruc- 
tures (|Chenl . 120051 : IZimenl . l200ll ). For compelling recent 
developments and advances in this latter are a we refer 
the in teres ted readers to the rece nt surveys by iBalandinl 
([20051 ) and lBalandin et all (|2007l ). together with the orig- 
inal literature cited therein. 

In this spirit, a primary building block for doing 
phononics is a setup that rectifies heat flow; i.e. a ther- 
mal rectifier/diode. Such a device acts as a thermal con- 
ductor if a positive thermal bias is applied while in the 
opposite case of a negative thermal bias it undergoes poor 
thermal conduction, thus effectively acting as a thermal 
insulator, or possibly also vice versa. The concept of such 
a thermal diode is sketched in Fig. [T] 

Thermal diode 




FIG. 1 (color). Sketch of the modus operandi of a thermal 
diode. When the left end of the diode is at a higher temper- 
ature as compared to its right counterpart heat is allowed to 
flow almost freely. In contrast, when the right end is made 
hotter in reference to the left end the transduction of heat 
becomes strongly diminished. 

The concept of a thermal diode involves, just as in 
its electronic counterpart, the presence of a symmetry 
breaking mechanism. This symmetry breaking is most 
conveniently realized by merging two materials exhibit- 
ing different heat transport characteristics. Historically, 
IStarrl (|1935l) working at Rensselaer Polytechnic Institute 
in New York built a junction composed of a metallic cop- 
per part which he joined with its cuprous oxide phase; 
thus proving the working principle of rectifying heat in 
such a structure. Starr's thermal rectifier is physically 
based on an asymmetric electron-phonon interaction oc- 



curring in the interface of the two dissimilar materials. 
There exist a plenty of such macroscopic rectifiers which 
function via the difference of the material response due to 
temperature bias and/or o ther externally applied cont rol 
fields such as strain, etc., ( Roberts and Walkerl l201lh . 

The focus in this colloquium will be on a thermal rec- 
tification scenario that is induced by phonon transport 
occurring on the nanoscale. The concept of such a ther- 
mal re ctifier for heat was put forward bv lTerraneo et al\ 
lj2002h . The authors therein proposed to use a three- 
segment structure composed of different nonlinear lat- 
tice segments. The underlying physical mechanism re- 
lies on the resonance phenomenon for the temperature 
dependent power spectrum vs. frequency as a result 
of the nonlinear lattice dynamics. Subsequent ly, it has 
been shown t hat a modified two-segment setup ( Li et all , 
l2005al . 12004) yields considerably improved rectification 
charac teristics as compared to the original three-segment 
setup (jTerraneo et "^I l2(K)2h . These pioneering works in 
turn ignited a flurry of activities, manifesting different 
advantageous features and characteristics. The theoreti- 
cal and numerical efforts culminated in a first experimen- 
tal validation of such a thermal rectifier in 2006: The de- 
vice itself is based o n an asymmetric nanotube structure 
(|Chang et aZl 120061) . The concept of this latter thermal 
diode together with its explicit experimental setup is de- 
picted with Fig. [2J 

The thermal diode presents an important first step to- 
wards phononics. For performing logic operations and 
useful circuitry, however, additional control mechanisms 
for heat are required. This comprises the task of devis- 
ing (i) the thermal analog of an electronic transistor, (ii) 
thermal logic gates and, as well, (iii) a thermal memory. 
The physical concept of these salient phononic building- 
blocks will be elucidated in Sect. II. The main physical 
feature that enters the function of such phonon devices is 
the occurrence of negative differential thermal resistance 
(NDTR); the latter being a direct consequence of the 
inherently acting nonlinear dynamics with its intriguing 
nonlinear response to an externally applied thermal bias. 

With Sect. Ill we investigate how to put these ther- 
mal phononic concepts to "action" by using realistic 
nanoscale structures. The actual operations of such 
devices rest on extended molecular dynamics simula- 
tions which serve as a guide for implementing its ex- 
perimental realization. The control of heat flow in the 
above mentioned phononic building blocks is managed 
mainly by applying a static thermal bias. More in- 
triguing control of transport emerges when the manip- 
ulations are made explicitly time-dependent or by use 
of different external forces such as the application of 
magnetic fields. As detailed with Sec. IV such manip- 
ulation scenarios then generate new roadways towards 
fine-tuned control and counterintuitive response behav- 
iors. Originally, such dynamic control has been imple- 
mented for anomalous particle transport by taking the 
system dynamics out-of-equilibrium: Doing so results in 
intriguing phenomena such as Brownian motor (ratchet- 
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FIG. 2 (color online). Concept and experimental setup of a 
nanoscale thermal diode. An asymmetric nanotube is placed 
over two electrodes with one serving as a heater. For fur- 
ther details on the experimental procedures and the exper- 
imental findings w e refer the readers to the experiment by 
IChang et all (2006) and the detailed discussion in Sect. Ill 
below. Adapted from lChang et all (120061 ). 



like) transport, absolu t e neg a tive mobility and alike 
(lAstumian and Hanggil 2002t lHanggi and Marchesonil 
120091 : lHanggi et all 1200511 . A similar reasoning can be 
put to work for shuttling heat in appropriately designed 
lattice structures, as detailed in Sect. IV. In Sect. V, we 
summarize our main findings, discuss yet some additional 
elements for phononic concepts and reflect on future po- 
tential and visions to advance the field of phononics from 
its present infancy towards a mature level. 



II. PHONONICS DEVICES: THEORETICAL CONCEPTS 
A. Thermal diode: Rectification of heat flow 

The task of directing heat for information processing 
as in electronics requires a toolkit with suitable building 
blocks, namely those nonlinear components that mimic 
the role of diodes, transistors, and alike, known from elec- 
tronic circuitry. A first challenge then is to design the 
blueprints for such components that function for heat 
control analogous to the building blocks for electronics. 



This objective is best approached by making use of the 
nonlinear dynamics present in anharmonic lattice struc- 
tures in combination with the implementation of a sys- 
tem inherent symmetry breaking. We start with the dis- 
cussion of the theoretical design for thermal diodes that 
rectify heat flow. 



1. Two-segment thermal diode 

In order to achieve thermal rectification, we exploit 
the nonlinear response mechanism as it derives from in- 
herent temperature-dependent power spectra. An every- 
day analog of such a nonlinear frequency response is a 
playground swing when driven into its large amplitude 
regime via parametric resonance. The response is opti- 
mized whenever the natural frequencies match those of 
the perturbations. Likewise, energy can be transported 
across two different segments when the corresponding vi- 
brational frequency response characteristics overlap. 

More precisely, whenever the power spectrum in one 
part of the device matches with its neighboring part, 
we find that heat energy is exchanged efficiently. In 
the absence of such overlapping spectral properties, the 
exchange of energy becomes strongly diminished. Par- 
ticularly, the response behavior of realistic materials is 
typically anharmonic by nature. As a consequence, the 
corresponding power spectra become strongly dependent 
on temperature, see in Appendix A, Subsect. 3. If an 
asymmetric system is composed of different parts with 
differing physical parameters the resulting temperature- 
dependence of the power spectra will differ likewise. 

Based on these insights, a possible working principle 
of a thermal diode goes as follows: If a temperature bias 
makes the spectral features of different parts overlap with 
each other, we obtain a favorable energy exchange. In 
contrast, if for the opposite temperature bias the spec- 
tral properties of the different parts fail to overlap ap- 
preciably a strong suppression of heat transfer occurs. 
In summary, this match/mismatch of spectral properties 
provides the salient mechanism for thermal rectification, 
see in Fig- El^b) and (c). 

Because the power spectra of an arbitrary nonlinear 
material typically become temperature-dependent, the 
use of any asymmetric nonlinear system is expected to 
display an inequivalent heat transport upon reversal of 
the temperature bias. It is, however, not a simple task 
to design a device that results in physical designated and 
technologically feasible thermal diode properties. After 
having investigated a series of possible setups we designed 
a thermal diode model that perfo rms effic i ently over a 
wide range of system parameters ( Li et al\ , [2004). The 
blueprint of this device consists of two nonlinear segments 
which are weakly coupled by a linear spring with strength 
feint . Each segment is composed of a chain of particles in 
which each individual particle is coupled with its nearest 
neighbors by linear springs. This whole nonlinear two- 
segment chain is each subjected to a cosinusoidal varying 
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FIG. 3 (color online). Concept of a thermal diode. Panel (a). 
Blueprint of an efficient two-segment thermal diode composed 
of two different Frenkel-Kontorova chains. The left (red) seg- 
ment consisting of a chain of particles subjected to a strong 
cosinusoidal varying on-site potential (illustrated by the large 
wavy curve) is connected to the right (blue) segment pos- 
sessing a relatively weak on-site potential, (b) In the case 
that the temperature XL in the left segment is colder than 
the corresponding right temperature Tr, i.e. Tl < Tr, the 
power spectrum of the particle motions of the left segment is 
weighted at high frequencies. This is so because of the dif- 
ficulty experienced by the dwelling particles to overcome the 
large barriers of the on-site potential. In contrast, the power 
spectrum of the right segment is weighted at low frequencies. 
As a result, the overlap of the spectra is weak, implying that 
the heat current J becomes strongly diminished, (c) Here the 
situation is opposite to panel (b). With Tl > Tr the parti- 
cles can now move almost freely between neighboring barriers. 
Consequently the power spectrum extends to much lower fre- 
quencies, yielding an appreciable overlap with the right placed 
segment. This in turn causes a sizable heat current. Panel 
(d). Heat current J vs. the relative temperature bias A, 
as defined in the inset, for three diffe rent values of the ref - 
ere nce temperature To. Adapted from IWang and Tl (2008a) 
and lLi et al\ (|2004f ). 
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on-site potential; the latter is provided by the coupling 
to a substrate. These individual chains are therefore de- 
scribed by a Frenkel-Kontorova (FK) lattice dynamics, 
cf. Appendix. The scheme of this thermal diode is de- 
picted in Fig. [3ja). 

The key feature of this FK-diode setup is the chosen 
difference in the strength of the corresponding on-site 
potential. At low temperature, the particles are confined 
in the valleys of the on-site potential. Thus the power 
spectrum is weighted in the high frequency regime. At 
high temperatures the particles assume su fficiently large 
kinet ic energies so that thermal activation (jHanggi et all 
1990) across the inhibiting barriers becomes feasible. The 
corresponding power spectrum is then moved towards 
lower frequencies. By setting the strength of the on-site 
potential in the two segments at distinct different lev- 
els, see Fig. we achieve the desired strong thermal 
rectification. Note that the barrier height of the on-site 
potential for the right segment is chosen sufficiently small 
so that the corresponding particles are allowed to move 
almost freely, both in the low and in the high tempera- 
ture regime. In the case that the left end is set at the low 
temperature, its power spectrum is weighted within the 
high frequency regime. This in turn causes an apprecia- 
ble mismatch with the right segment, see Fig.[3Jb). In the 
opposite case, when the left end is set at the high temper- 
ature, its weighted power spectrum moves towards lower 
frequencies, thus matching considerably with the right 
segment, see Fig. IHJc). 

In Fig. Eld) the resulting stationary heat current J 
(see in Appendix A) vs. the relative temperature bias 
A is depict ed for three va lues of the reference temper- 
atures To (|Li et all 12004 ) . The relation between the 
dimensionless temperature and the actual physical tem- 
perature can be found in the Appendix as well. It is 
shown that when A > (i.e. Tl > Tr), the heat current 
gradually increases with increasing A, i.e. the setup be- 
haves as a 'good' thermal conductor; in contrast, when 
A < (Tl < Tr), the heat current remains small. The 
two-segment structure thus behaves as a 'poor' thermal 
conductor, i.e. it mimics a thermal insulator. 

For a given setup, the heat current through the system 
is mainly controlled by its interface coupling strength 
fei n t . Figure f?] depicts the temperature profiles for dif- 
ferent fei n t and for two oppositely chosen bias strengths. 
There occurs a large temperature jump at the interface. 
The size of the jump is larger for negative bias A (filled 
symbols in Fig. @| than for positive bias A. In the case 
with negative bias the temperature gradient inside each 
lattice segment almost vanishes; implying that the result- 
ing heat current is very small. This behavior is opposite 
to the case with positive bias (open symbols in Fig. 2]). 



2. Asymmetric Kapitza resistance 

The interface thermal resistance (ITR), also known 
as the Kapitza resistance, measures the interfacial re- 
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FIG. 4 (color online). Temperature profiles within a FK-FK 
thermal diode. The temperature profiles for various interface 
elastic constant ki nt = 0.01 (square), 0.05 (circle) and 0.2 (tri- 
angle). To = 0.07, N = 100. The open symbols correspond 
to a positive temperature bias A = 0.5, while the filled sym- 
bols correspond to negative temperature bias A = —0.5. The 
dashed line indicates the interface. The temperature jumps 
at the interface are due to the Kapit za resistance w hich is 
addressed in Sec. II. A. 2. Adapted from|L i et al\ (|2004l ). 



sistance to heat flow dPollackl I1969L ISwartz and Poll 
Il989f ). It is defined as, R = AT/ J, where J is the 
heat flow (per area) and AT denotes the temperature 
jump between two sides of the interface. The origin of 
the Kapitza resistance can be traced back to the hetero- 
geneous electronic and/or vibrational properties of the 
different materials making up the interface at which the 
energy carriers become scattered. The amount of rela- 
tive transmission depends on the available energy states 
on either of the two sides of the interface. Originally, 
this phenomenon was discovered by Kapitza i n 1941 in 
exper iments detecting super-fluidity of He II ( Kapitzal 
I1941D . 

The high thermal rectification in the above discussed 
model setups is mainly due to this inter face effec t . In or - 
der to further improve the performance iLi et all (2005a) 
studied the ITR in a lattice consisting of two weakly 
coupled, dissimilar anharmonic segments, exemplified by 
a (FK)-chain segment and a neighboring Fermi-Pasta- 
Ulam (FPU)-chain segment. 

Not surprisingly, the ITR in such a setup depends on 
the direction of the applied temperature bias. A quan- 
tity that measures the degree of overlap of power spectra 
between left(L) and right(R) segments reads 



5' 



f °° P L {w)P R {u>)du 



(1) 



numerical simulations then reveal that 
(S+/S-) Sr with S R = 1.68 ± 0.08, and 

\J+/J-\ ~ (S+/S-) 6 '', with Sj = 1.62 ±0.10. The nota- 
tion +/— indicates the case with a positive thermal bias, 
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A > 0, and a negative thermal bias, A < 0, respectively. 
These findings thus support the strong dependence of 
thermal resistance on this very overlap of power spectra. 

The physical mechanism of the asymmetric ITR be- 
tween dissimilar anharmonic lattices can therefore be re- 
lated to the match/mismatch of the corresponding power 
spectra. As temperature increases, the power spectrum 
of the FK lattice shifts downwards, towards lower fre- 
quencies. In contrast, the power spectrum of the FPU- 
segment, however, shifts upwards, towards higher fre- 
quencies. Due to these opposing shifts, it becomes ev- 
ident that upon a reversing the thermal bias the amount 
of match/mismatch in such a FK-FPU setup can even 
exceed that of the above considered FK-FK setup. Con- 
ceptually this results in an even stronger thermal rectifi- 
cation. 



B. Negative differential thermal resistance: The thermal 
transistor 

The design and the experimental realization of the 
thermal diode presents a striking first step towards the 
goal of doing phononics. A next challenge to be overcome 
is then a design for a thermal transistor, which allows for 
an a priori control of heat flow much alike the familiar 
control of charge flow in a Field Effect Transistor. Like 
its electronic counterpart, a thermal transistor consists of 
three terminals: The drain (D), the source (S), and the 
gate (G). When a constant temperature bias is applied 
between the drain and the source, the thermal current 
flowing between source and drain can be fine-tuned by 
the temperature that is applied to the gate. Most impor- 
tantly, because the transistor is able to amplify a signal, 
the changes in the heat current through the gate can in- 
duce an even larger current change from the drain to the 
source. 

Towards the eventual realization of such a thermal 
transistor device we consider the concept depicted with 
Fig. EJa). It uses a one-dimensional anharmonic lattice 
structure where the temperature at its ends are fixed at 
Td and Tg (Tb > Is), respectively. An additional, third 
heat bath at temperature To controls the temperature 
at node O, so as to control the two heat currents Jd and 
Js = Jd + Jo- Let us next define the current amplifica- 
tion factor a. This quantity describes the amplification 
ability of the thermal transistor, as the change of the heat 
current Jd (or Js, respectively), divided by the change 
in gate current Jo , which acts as the control signal. This 
amplification quantity explicitly reads 



to yield 



ks/(js + r D )| 



(5) 



a = |<9J D /<9Jo| 



(2) 



Eq. @ can be recast in terms of the differential thermal 
resistance of the segment S, i.e., 

r s = (dJ s /dT )-' =const (3) 

and that of the neighboring segment D, i.e., 

r D = - (dJv/dT ) T ^ const (4) 



It can readily be deduced that if rg and ro are both pos- 
itive then a emerges to be less than unity. Consequently, 
such a thermal transistor cannot work. 




FIG. 5 (color online). Concept of a thermal transistor. Panel 
(a). A one-dimensional lattice is coupled at its two ends to 
heat baths at temperature Td and Ts with Td > Ts. A third 
heat bath with temperature Ts < To < Td can be used to 
control the temperature at the node O, so as to control the 
heat currents Jd and Js. Panel (b). In an extended regime, 
as the temperature To increases, the thermal bias (Td — To) 
at the interface decreases while the power spectrum in the left 
segment of the control node (wine-color) increasingly matches 
with the power spectrum of the neighboring segment that is 
connected to the drain D (red). This behavior is depicted 
with the insets where the corresponding power spectra of the 
left-sided and right-sided interface particles are depicted for 
three representative values of To, as shown by the three ar- 
rows. The resulting drain current Jd increases with increas- 
ing overlap until it reaches a maximum, and starts to decrease 
again. 

For a transistor to work it is thus necessary that the 
current amplification factor obeys a > 1. This implies 
a negative differential thermal resistance (NDTR); i.e., 
it requires a transport regime wherein the heat current 
decreases with increasing thermal bias. Such behavior 
occurs with the thermal diode characteristics depicted 
with Fig. [3] (d) ; note the behavior for the case with filled 
(blue) triangles with A increasing from ~ [—0.5,-0.2]. 
It should be pointed out that such an (NDTR)-behavior 
is in no conflict with any physical laws 1 . 

While negative differential electric r esistance has long 
been realized and extensively studied (Esaki, 1958J), the 



1 Note that this NDTR should not be confused with absolute neg 
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concept of NDTR has been proposed more recently only 
(|Li et aLl . [200ll2006al) . With temperature T Q increasing, 
the sensitively temperature-dependent power spectra of 
the two segments match increasingly better which not 
only offsets the effect of a decreasing thermal bias {To — 
To) but even induces an increasing heat current. This 
behavior is illustrated in Fig. \5[b). 

A system displaying NDTR constitutes the main ingre- 
dient for operation of a thermal transistor. The scheme 
of a thermal transistor is shown in Fig. |6ja). In order to 
make this setup physically more realistic, a third segment 
(G) is connected to the node O. This is done so, because 
in an actual device it is difficult to directly control the 
temperature of node O, being located inside the device. 
Using different sets of parameters, this thermal transistor 
can work either as a thermal switch, Fig.[6jb), or also as 
a thermal modulator, Fig. [IJc). 

The key prerequisite for a thermal transistor, i.e. the 
NDTR phenomenon has been investigated in various 
other system s recently , e.g.. for high dimensional lat- 
tice models ( Lo et all . 120081 ). A gas- liquid transition 
has als o been utilized in the d esign of a thermal tran- 
sistor ( Komatsu and Itol . 1201 ll ). The condition for the 
existence of a NDTR regime is more stringent than 
that of obtaining merely thermal rectification. The 
crossover from existence to nonexistence of NDTR have 
been investigated for a set of different la ttice structures , 
both by an a lytica l and numerical means (iHe et aLl . l2009t 
IShao et all l2009h . Because the NDTR in these lattice 
models basically derives from an interface effect, it is 
plausible that with a too large interface coupling strength 
feint or a too long lattice length, it is rather the ther- 
mal resistance of the involved segments than the inter- 
face resistance that rules the NDTR. As a consequence, 
the NDTR effect typically becomes considerably sup- 
pressed. We therefore expect that NDTR can experi- 
mentally be realized with nano-scale materi als, for exam- 
ple u s ing nano-tu bes, nano- wires, etc., see ( Chang et all 
2007. l2008L[2006h . This issue is addressed in greater de- 
tail with Sect. Ill below. 



C. Thermal logic gates 

The phenomenon of NDTR provides not only the func- 
tion for thermal switching and thermal modulating, but 
also the essential input in devising thermal logic gates. 
Setups for all major thermal logic gates able to per- 



ative thermal resistance around an equilibrium working point at 
thermal bias A = 0. The latter would imply that heat flows 
from c old to hot wh ich violates the principle of Le Chatelier- 
Brown llCallenl . H960h . i.e. no opposite response to a small exter- 
nal perturbation around a thermal equilibrium is possible. Such 
an anomalous response behavior is, however, feasible when the 
syste m is taken (at a zero bias) into a stationary non-equilibrium 
state llHanggi and Marchcsoni, 2 0091) . 
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FIG. 6 (color online). Thermal transistor, (a) Sketch of a 
thermal transistor device. Just as in the case of an electronic 
transistor, it consists of two segments (the Source and the 
Drain) and, as well, a third segment (the Gate) where the 
control signal is injected. The temperatures To and Ts are 
fixed at high, T+, and low, T_, values. The negative dif- 
ferential thermal resistance (NDTR) occurs at the interface 
between O and O'. This part then makes it possible that over 
a wide regime of parameters, when the gate temperature Tg 
rises, not only Js but also Jd may increase. Using different 
system parameters allow then different function, this being 
either a thermal switch (panel (b)) or also a thermal modu- 
lator (panel (c)). Panel (b). Function of a thermal switch: 
At three working points where Tg = 0.04, 0.09 and 0.14, the 
heat current Jd equals Js, yielding Jg = Js — Jd vanishing 
identically. These three working points correspond to (stable) 
"off", (unstable) "semi-on" and (stable) "on" states, at which 
Jd differs substantially. We can switch, i.e., forbid or allow 
heat flowing by setting Tg at these different values. Panel 
(c). Function of a thermal modulator. Over a wide tempera- 
ture interval of gate temperature Tg, depicted via the hatched 
working region, the heat current Jg remains very small, i.e. 
it remains inside the hatched regime, while the two heat cur- 
rents Js and Jd can be co ntinuous l y contr olled from low to 
high values. Adapted from lLi et al\ |2006a). 
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form logic operation s have been put forward recently 
(IWang and Lil . l2007l) . 

In digital electric circuits, two boolean states "1" and 
"0" are encoded by two different values of voltage, while 
in a digital thermal circuit these boolean states can be 
defined by two different values of temperature T on and 
T Q ff. In the following we discuss the way to realize such 
individual logic gates. 

A most fundamental logic gate is the signal repeater 
which 'digitizes' the input. Namely, when the input tem- 
perature is lower/higher than a critical value T c , with 
T ff < T c < T on , the output is set at T Q g/T on , respec- 
tively. This does not present a simple task; it must take 
into account that small errors may accumulate, thus lead- 
ing eventually to incorrect outputs. 

The thermal repeater can be obtained by use of ther- 
mal switches. Let us inspect again the thermal switch 
shown in Fig. |5ta), in which the Tc-dependence of the 
heat currents Jd , Js and Jq is illustrated with Fig. [SJb) . 
When the gate temperature Tq is set very close but not 
precisely at T g/T on , then the heat current in the gate 
segment makes the temperature in the junction node O 
approach more closely T g/T on . Therefore, upon con- 
necting such switches in series, which involves plugging 
the output (from node O) of one switch into the gate 
of the next one, the final output will asymptotically ap- 
proach that of an ideal repeater, i.e. T ^/T on , whichever 
is closer to the input temperature Tq. 

A NOT gate reverses the input; it yields the response 
"1" whenever it receives "0", and vice versa. This re- 
quires that the output temperature falls when the input 
temperature rises, and vice versa. This feature can be 
realized by injecting the signal from the node G and col- 
lecting the output from the node O', cf. Fig. [6ja). The 
NDTR between the nodes O and O' again plays the key 
role. A higher temperature Tq induces a larger thermal 
current Jd and therefore increases the temperature bias 
in segment D. To- thus decreases (notice that To is fixed) 
and a negative response is thus realized, cf. Fig. EJa). 
Suppose To- equals Tgf and Tg?, with T°?>T ?, when 
Tq equals T Q g and T on , respectively. A remaining prob- 
lem is that both Tg? and Tg? are higher than T c (in fact 
even higher than T on ) . This in turn will be always treated 
as a logical "1" by the following device. This problem can 
be solved if we apply a "temperature divider" , the coun- 
terpart of a voltage divider, which is depicted in inset 
of Fig. [7] (b). Its output equals a ratio of its input. By 
adjusting this ratio, we can make the output of the tem- 
perature divider be higher/lower than T c when its input 
equals Tq?/T™. We then employ a thermal repeater to 
digitize the output from the temperature divider. In do- 
ing so the function of a NOT gate is therefore realized, 
as depicted with Fig. [7] (b). 

Similarly, an AND gate is a three-terminal (two inputs 
and one output) device. The output is "0" if either of the 
inputs are "0" . Because we now have the thermal signal 
repeater at our disposal, an AND gate is readily realized 
by plugging two inputs into the same repeater. It is clear 
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FIG. 7 (color online). Negative response and thermal NOT 
gate, (a) Temperature of the node O' as a function of Tq 
for the setup shown in Fig. [6ja). In a wide range of Tg, To' 
decreases as Tq increases, (b) Function of the thermal NOT 
gate. The thin (blue) line indicates the function of an ideal 
NOT gate. Inset: Structure of a two-resistor voltage divider, 
the counterpart of a temperature divider, which supplies a 
voltage lower than that of the battery. The output of the 
voltage divid e r is: V ou t = VR2/{R\ + R2)- Adapted from 
IWang and Ll (|2007l ). 



that when both inputs are "I" , then the output is also 
"I" ; and when both inputs are "0" , then the output is 
also "0" . By adjusting some parameters of the repeater, 
we are able to make the final output to be "0" when 
the two inputs are opposite, therefore realizing a thermal 
AND gate. An OR gate, which exports "I" whenever 
the two inputs are opposite, can also be realized in the 
similar way. 



D. Thermal memory 

Towards the goal of an all phononic computing, yet an- 
other indispensable element, besides thermal logic gates, 
are thermal memory elements that enable the storage of 
information via its encoding by heat or temperature. 

A possible such setup acting a s a therm al memory has 
been proposed by IWang and Lil ( 2008br l. Its blueprint 
has much in common with the working scheme for the 
thermal transistor, cf. Fig. [^a). Adjusting param- 
eters appropriately, negative differential thermal resis- 
tance can be induced between the chain segments con- 
necting to node O and node O', respectively. With fixed 
Ts and Tb, obeying Tg < 2d, and a heat bath at tem- 
perature Tq that is coupled to the node G, three possible 
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working points for Tq, i.e, T q, T sem i_ on and T on can be 
realized. At these three working points the gate cur- 
rent vanishes, Jq — 0, thus balancing perfectly Js and 
Jd- Upon analyzing the slope of Jg at these very points 
two of these, i.e., (T g,T on ) denote stable working points 
and the intermediate one, T som i-oii, is unstable. Notably, 
these working states remain stationary when the heat 
bath that is coupled to terminal G is removed. Now, 
however, the corresponding temperatures at these work- 
ing points exhibit flu ctuations. As is we ll known from 
stochastic bistability (jHanggi et all , Il990|) , small fluctu- 
ations around these working points drive the system con- 
sistently back to the stable working points and away from 
the unstable working point. Accordingly, the system pos- 
sesses two long lived meta-stable states, To = T on and 
To=T oS , while T = T somi _ on is unstable. 

The stability of these so adjusted thermal states at T on , 
T fj implies that these states remain basically unchanged 
over an extended time span in spite of thermal fluctua- 
tions. This holds true even if a small external pertur- 
bation, as for example imposed by a small thermometer, 
reading the temperature at the node O, is applied. 

The working principle of a Write-and-Read cycle of 
this so designed thermal memory is depicted with Fig. [8] 
Starting out at time t = from a random initial prepa- 
ration of all particles making up the memory device the 
local temperature at site O relaxes to its stationary value 
To ~ 0.18 (initializing stage). For the case in panel (b) 
in Fig. [5] the Writer, being prepared at its boolean tem- 
perature value T g is next connected to the site O (writ- 
ing stage) . As can be deduced from panel Fig. [HJb) , the 
temperature To quickly relaxes during this writing cy- 
cle to this boolean value, and maintains this value over 
an extended time span, even after the Writer is removed 
(maintaining stage). More importantly, To self-recovers 
to this setting temperature T g after being exposed to the 
small perturbation as induced by the Reader (a small 
thermometer) during the following data-reading stage. 
The data stored in the thermal memory are therefore 
precisely read out without destroying the memory state. 

In panel Fig. [5(c) , the Write-and-Read process corre- 
sponding to the opposite boolean writing temperature 
value, i.e. To = T on , is depicted. Accordingly, this so 
engineered Write-and-Read cycle with its two possible 
writing temperatures T g and T on does accomplish the 
task of a thermal memory device. 



III. PUTTING PHONONS TO WORK 

In the foregoing Sec. II, we have investigated vari- 
ous setups involving stylized nonlinear lattice structures 
in order to manipulate heat flow. We next discuss how 
to put these concepts into practical use with physically 
realistic systems. In doing so we consider numerical stud- 
ies of suitably designed nanostructures which exhibit the 
designated thermal rectification properties. This discus- 
sion is then followed with first experimental realizations 
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FIG. 8 (color). Write-and-Read cycle in a thermal memory, 
(a) Illustration of the four stages of a complete writing-reading 
process. The left and right ends of this thermal memory are 
coupled to heat baths at fixed temperatures T = 0.03 and 
T — 0.3, respectively. During the initializing stage the sys- 
tem relaxes towards its stationary state with the tempera- 
ture at site O being To = 0.18. The different stages during 
this cycle are indicated by separating dashed vertical lines 
in panels (b) and (c). (b) Time evolution of the local tem- 
perature To(t) (dotted wine color) of the central particle at 
site O, the temperature of the Writer T wr iter(£) (solid olive) 
and the reader T rea d e r(£) (dashed violet) during a complete 
Write-and-Read process. The time-dependent local tempera- 
ture of the particle (with unit mass) at site O is obtained by 
integrating and averaging the squared velocity vo(t) over a 
short time span (10 4 time steps), centered at time-instant t. 
Writer and Reader are realized with short lattice segments, 
while T wr iter and T rea der are the averaged temperatures of all 
particles within the Writer and Reader, respectively. Dur- 
ing the writing stage the node O is coupled to one end of the 
Writer whose opposite end is coupled to a heat bath set at the 
"off"-value T ff = 0.04. During the reading stage, the reader 
set at initial temperature 0.11 is attached to the node O. (c) 
The opposite Write-and-Read cycle with the Writer set at 
the "on " -value with T on = 0.27. Adapted from IWang and Lil 
(|2008b| y 



of a thermal diode and the thermal memory. 

Among the plenty of physical materials that come 
to mind, low dimensional nanostructures like nan- 
otubes, nanowires and graphene likely offer optimal 
choices to realize the desired thermal rectification fea- 
tures obtained with nonlinear lattice studies. It is 
known from theoretical studies (|Dharl . 120081 : iLepri et all , 
120031: ISaito and Dharl.l2010l) and experi mental validation 
(jChang et all l2008t iGhosh et all , I2010T ) , that the charac- 
teristi cs of heat flow (such as the valid i ty of the Fourier 
Law (|Dharl . 120081: ILepri et all , 120031: ISaito and Dhal 
l2010h ) sensitively depends on spatial dimensionality and 
the absence or presence of (momentum) conservation 
laws. 
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This is even more the case for nano materials, wherein 
due to the l imite d size, the discrete phonon spectrum 
(|Yang et "flll. l2010h results in a distinct dependence of the 
thermal quantities on the specific geometrical configura- 
tion, mass distribution and ambient temperature. Over- 
all, this makes nanosized materials promising candidates 
for phononics. Because the experimental determination 
of thermal transport properties is not straightforward the 
combined use of theory and numerical simulation is indis- 
pensable in devising phononics devices. Moreover, novel 
atomistic computational algorithms have been developed 
which faci litate the study of experimentally relevant sys- 



wmcn iacuitatc tnc study < 
tern sizes ( Li et all |2010| ). 



A. Thermal diodes from asymmetric nanostructures 
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FIG. 9 (color online). Schematic system setups for ther- 
mal diodes devised from asymmetric nano structures: (a) 
a typical (n, 0)/(2n, 0) carbon nanotube junction; (b) a 
carbon-nanocone based thermal rectifier; (c) two asymmet- 
ric graphene nanoribbons (GNRs), a trapezia-shaped GNR 
and, as well, a structure made of two-rectangular GNRs of 
different width; (d) a topological configuration based on a 
Mobius graphene stripe. The red parts denote the corre- 
sponding heat bath regions. The bath temperature of two 
ends are denoted as Tt (Top) and Tb (Bottom) in (b) and (c), 
and as Tl (Left) and Tr (Rig ht) in (a) and (d), r espectively. 
These f igures are adapted f romlWu and Lil (l2007f ). lYang et al\ 
(l2008l h lYairg et al\ |2009r i, |jiang et al\ (|2010T) . 

Carbon nanotubes (CNTs) have recently attracted at- 
tention for applications in nanoscale electronics, mechan- 
ical and thermal devices. CNTs posses a high the rmal 
conductivity at room temperature ( Kim et a/. . l200ll) and 
phonon mean free paths that extend over the length scale 
of structural rip ples; thus providing ideal phonon waveg- 
uide properties ( Chang et all 120071) . 

Let us consider thermal rectification in s ingle- walled 
carbo n nanotube (SWCNT) based junctions (|Wu and Lil . 
120071 ). A typical (n, 0)/(2n, 0) intramolecular junction 
structure is shown in Fig. m which the structure con- 
tains two parts, namely, a segment with a (n, 0) SWCNT 
and a segment made of a (2n, 0) SWCNT. These two seg- 
ments are connected by m pairs of pentagon-heptagon 
defects. By use of non-equilibrium molecular dynamics 



(NEMD) simulations, one finds that the heat flux from 
the (2n, 0) to (n, 0) tube exceeds that from the (n, 0) 
tube to the (2n, 0) segment. The corresponding thermal 
rectification increases upon raising the temperature bias. 

Another beneficial feature is that the rectification is 
only weakly dependent on the detailed structure of the 
interface, assuming that the connecting region is suffi- 
ciently short. This is so because heat is predominantly 
carried by long wavelength phonons, which are scattered 
mainly by large defects. Just as with nonlinear lat- 
tice models, the match/mismatch of the energy spec- 
tra around the interface is the underlying mechanism 
for rectificat i on. F urthermore, in the elongated structure 
(IWu and Lil . |2008) the heat flux becomes smaller than 
that of the non-deformed structure; its temperature de- 
pendence, however, becomes more pronounced, due to 
the intrinsic tensile stress. 

The carbon nanocone depicted in Fig. [HJb) is yet an- 
other carbon based material exhibiting a high asymmet- 
ric geometry. Its thermal properties were investigated by 
I Yang et~al\ ( 20081 ). A temperature difference, A, between 
the two ends of nanocone is introduced, being positive 
when the bottom of nanocone is at a higher tempera- 
ture, and is negative in the opposite case. It was found 
that the nanocone behaves as a "good" thermal conduc- 
tor under positive thermal bias and as a "poor" thermal 
conductor when exposed to a negative thermal bias. This 
suggests that the heat flux runs preferentially along the 
direction of decreasing diameter. 

In order to compare the impact of mass-asymmetry 
versus geometry-asymmetry for thermal rectification a 
nanocone structure with a g raded mass distribution was 
discussed ( Yang et all [2008) . The mass of the atoms of 
the nanocone are devised to change linearly; that is, the 
top atoms possess a minimum mass Men and the bot- 
tom atoms are at the maximum mass of AMq 12, with 
Mci2 being the mass of 12 C atom. NEMD results then 
yield a rectification ratio, | J+ — J-|/| J-\, that is 10% for 
a nanocone with uniform masses and is 12% for the mass- 
graded nanocone (at identical thermal bias A). This 
2% increase for the mass-graded distribution evidences 
that the role of geometric asymmetry is more effective in 
boosting thermal rectification. 

These geometric asymmetric SWCNT junctions all ex- 
hibit thermal rectification. Compared to the SWCNT 
setup, it is easier, however, to control the shape of 
graphene by the nano cutting technology such as he- 
lium ion microscope. Graphene nanoribbons (GNRs), cf. 
Fig. [9jc), present promising element s for nanoel e ctron - 
ics and nanophononics. For instance, I Yang et~al\ (2009) 
have demonstrated tunable thermal conduction in GNRs. 
They studied the direction dependent heat flux in asym- 
metric structural GNRs, and explored the impacts of 
both, GNR shape and size, for the rectification ratio. 
Two types of GNRs were considered: trapezia-shaped 
GNR and two-rectangular GNRs of different width, Fig. 
He). 

These two types of GNRs behave as a "good" thermal 
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conductors under positive thermal bias (i.e. bottom end 
at higher temperature) and as a "poor" thermal conduc- 
tor under negative thermal bias (i.e. top end at higher 
temperature). This finding is similar to the rectifica- 
tion phenomena o bserved in carbon nanocone structures 
(|Yang et adl2008h . which again results from the amount 
of match/mismatch between the respective phonon spec- 
tra. 

It is interesting to note that the rectification ratios of 
GNRs are substantially larger than those in nanocone 
and SWCNT junctions. Wi th A=0.5, the re ctification 
ratio of the nanocone is 96% ( Yang et a/.l . l2008h and that 
of SWCNT intram olecular junction is only around 15% 
(TWu and Lil . l2007h . while the rectification ratio in GNRs 
is about 270% and 350% for two-rectangular GNRs and 
the trapezia-shaped GNR, respectively. The rectifica- 
tion ratio of trapezia-shaped GNR is thus larger than 
that of two- rectangular GNR under the same thermal 
bias difference. This feature is consistent with the phe- 
nomenon that the carbon nanocone (i.e. the geomet- 
ric graded structure) possesses a larger rectification ra- 
tio than the two-segment carbon nanotube (n, 0)/(2n, 
0) intra-molecular junction. The phonon spectra change 
continuously in the graded structures which implies a 
more efficient control of heat flux. Similar rectifica- 
tion features we re observed in graphene nanoribbons by 
iHu et al\ (I2009D . 

Yet another advantage of a GNR based thermal rec- 
tifier is its weak dependence on length. Because en- 
ergy is transported ballistically in graphene, the heat 
flux is essentially independent on size. Both J + and 
J_ are saturated when the GNR length is longer than 
~ 5.1 nm, yielding that the rectification ratio remains 
constant around 92%. Moreover, compared to a setup 
composed of a single-layer graphene, an even larger rec- 
tification ratio can be achieved in few-layer asymmetric 
structur es, this being a conseque nce of layer-layer inter- 
actions ( Zhang and Zhang] . l201l[ ). 

Thermal rectification can also be realized with a topo- 
logical graphene se tup such as the Mobius graphene strip 
([Jiang et all 120101) . Fig. E^d). The advantage of this 
topology induced thermal rectification is that it is prac- 
tically insensitive to temperature and size of the system. 
In this structure, the asymmetry originates from the in- 
trinsic topological configuration. 

Because phonons are strongly scattered due to the mis- 
match in vibrational properties of the materials forming 
the interface, which in addition depends on the sign of ap- 
plied thermal bias, this produces an asymmetric Kapitza 
contact resistance. Capitalizing on this idea a silicon- 
amorphou s poly ethylene thermal rectifier was designed 
(IHu et all I2008h . The result is that the heat current 
from the polymer to the silicon is larger than vice versa. 
To examine the origin of the thermal rectification, the 
density of states of phonons on each side of the interface 
was investigated. The phonon density of states signifi- 
cantly softens in the ploymer as it becomes warmer. This 
increases the density of states in the polymer at low fre- 



quencies. Those low frequency acoustic modes in silicon 
thus increase their transmission probability, yielding an 
enhanced thermal transport. 

The setups discussed thus far are all operating in 
steady state. A transient thermal transport exhibiting 
a time-dependent thermal rectifica tion has been investi- 
gated for a Y-SWCNT junction by iGonzalez Nova et al\ 
( 20091 ). From their molecular dynamics simulation it is 
reported that a heat pulse propagates unimpeded from 
the stem to the branches, undergoing little reflection. For 
the reverse temperature bias, however, there is a substan- 
tial reflection back into the branches. The transmission 
coefficient from the stem to the branches is more than 
four times that of the reverse direction. 

Let us reflect again on the choice of materials suitable 
in obtaining those theoretically predicted thermal recti- 
fication features. A thermal rectifier constitutes a two- 
terminal thermal device whose working principle rests 
upon the different overlap of the temperature dependent 
phonon spectra. For this to occur strong anharmonicity 
plays an important factor. Given the above discussion 
of geometric asymmetric structures a strong intrinsic an- 
harmonicity corresponds to diffusive phonon transport. 
However, diffusive phonon transport presents not a nec- 
essary condition for thermal rectification. In the asym- 
metric GNRs, phonon transport is almost ballistic. The 
conclusion is that thermal rectification exists for both, 
ballistic and diffusive phonon transport, assuming that 
a match/mismatch of the phonon spectrum of the two 
ends is present. This difference in phonon spectra can 
originate from different sources, as for example from a 
asymmetric mass distribution, different geometry, size or 
spatial dimension. 

Despite an abu ndance o f parameter dependent rec- 
tification s tudie s (IHu et all [20061: I Wu and Segall 120091: 
IWu et all 120091: lYang et all 120071 ). these were rarely 
tested against experiments. All those NEMD-studies 
rely on mathematically idealized material Hamiltonians, 
which in practice may still be far from physical reality. 
Only in situ experiments thus present the ultimate test 
bed for a validation of the wealth of available theoretical 
results and, even more importantly, dictate the necessary 
next steps towards an implementation of phononics. 



B. In situ thermal diodes from mass-graded nanotubes: 
Experiment 

On the experimental side, a pioneer th e rmal diode 
work has been performed by IChang et all (|2006f) . In 
their experiment, cf. Fig. [Kjl carbon nanotubes (CNTs) 
and boron nitride nanotubes (BNNTs) were gradually de- 
posited on the surface with the heavy molecules located 
along the length of the nanotube in order to establish an 
asymmetric mass distribution. As an important test, it 
has been experimentally checked that in unmodified NTs 
with uniform mass distribution the thermal conductance 
is indeed independent of the direction of the applied ther- 
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FIG. 10 (color online). Experimental realization of thermal 
rectification in nanotube junctions, (a) Scanning-Electron- 
Microscope images of boron nitride nanotubes (BNNTs) after 
deposition of CgHiaPt. (b) Graphical representation of the 
temperature changes of the heater (ATh) an d sensor (AT S ) 
for the nanot ubes before and aft er deposition of CgHiePt. 
Adapted from lChang et all (120061 ). 



mal bias. In clear contrast, however, the inhomogeneous 
NT-system in fact does exhibit asymmetric axial ther- 
mal conductance, with greater heat flow in the direction 
of decreasing mass density. The thermal rectification ra- 
tios, i.e. \J+ — J-\/\J-\, were measured as 2% and 7% for 
CNT and BNNT, respectively. The higher thermal rec- 
tification found in BNNT might originate from stronger 
nonlinearity, as induced by the ionic nature of the B-N 
bonds, which in turn favors thermal rectification. 

In order to understand the mechanism of rectifi- 
cation in Chang et al's experiment, subsequent work 
( Yang et "all . |2007( ) has studied the thermal properties in 
a on e-dimensional anh armonic lattice with a mass gradi- 
ent. lYang et all ((2007) found that in the ID mass-graded 
chain, when the heavy-mass end is set at high tempera- 
ture, the heat flux is larger than that under the reverse 
temperature bia s. This is con s istent with Chang's exper- 
imental results ( Chang et all [20061) . Also, the larger is 
the mass gradient the larger is thermal rectification. As 
mentioned before, this can be explained from differences 
in overlap of the respective phonon power spectra. 



C. Solid-state-based thermal memory: Experiment 

Apart from the experimental demonstration of a 
nanoscale thermal rectifier, a thermal memory device has 
recently also brought into operation experimentally by 
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FIG. 11 (color online). Experimental realization of the ther- 
mal memory, (a) Scanning-Electron-Microscope image of the 
memory device, consisting of a VO2 nanobeam connecting 
the input terminal (left side) and output terminal (right side). 
An equivalent thermal circuit is also depicted in the image. 
Panel (b). Output temperature, T out , as a function of input 
temperature, Tj n , upon heating (filled symbols) and cooling 
(open symbols). The characteristics of the thermal memory 
is tuned under three different voltage biases, as indicated in 
the panel. Panel (c). Write and Read out process of logical 
"1" (= HIGH temperature state) and "0" (= LOW temper- 
ature state). Upper part, Ti n (i) for writing (red line): To 
write a HIGH (=1) temperature state to the output terminal, 
a heating pulse generated by increasing the heating current 
was applied to the input terminal to set Ti n to 435K. To 
switch the output terminal to a LOW temperature state, a 
cooling pulse was applied to the input terminal to set T n to 
300K by natural cooling of Ti n to the substrate temperature 
Tbasc- Lower part, T out (t) for reading (blue line): When Ti n 
was returned to 360K for reading, the HIGH (LOW) state 
of ~ 332.1±0.1 K (~ 317±0.8K) was retained and could be 
read out at the output terminal T ou t. Here the voltage bias 
is 0.04V. Data and Figure adapted from lXie et all (|201lT l . 



IXie et all (f2PHll In their work, see in Fig. Illf a). a 
single-crystalline VO2 nanobeam is used to store and re- 
tain thermal information with two temperature states as 
input T; n and output T out , which serve as the logical 
boolean units '!"(= HIGH) and "0"(=LOW), for writ- 
ing and reading. This has been achieved by exploiting 
a metal-insulator transition. A nonlinear hysteretic re- 
sponse in temperature was obtained in this way. A volt- 
age bias across the nanobeam has been applied to tune 
the characteristics of the thermal memory. One finds that 
the hysteresis loop becomes substantially enlarged, see 
in Fig. ITTT b). and is shifted towards lower temperatures 
with increasing voltage bias. Moreover, the difference in 
the output temperature T out between its "HIGH" and 
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"LOW" temperature states increases substantially with 
increasing voltage bias. To demonstrate the repeatabil- 
ity of the thermal memory, they have performed repeated 
Write(High)-Read-Write(Low)-Read cycles using heating 
and cooling pulses to the input terminal, see panel (c) 
in Fig. QT] for more details. Repeated cycling over 150 
times proves the reliable and repeatable performance of 
this thermal memory. 



IV. SHUTTLING HEAT AND BEYOND 

The function of the various thermal electronic analog 
devices discussed in Sect. II used a heat control mech- 
anism which is based on a static thermal bias. In or- 
der to obtain an even more flexible control of heat en- 
ergy, being comparable with the richness available known 
for electronics, one may design intriguing phononic de- 
vices which utilize temporal, ac gating modulations as 
well. Such forcing makes possible the realization of a 
plenitude of novel phenomena such as the heat ratchet 
effect, absolute negative heat conductance or the real- 
i zation of Brownian (h e at) m o tors, to name but a few 
(lAstumian and Hanggil. l2002t lHanggi and Marchesonil 
120091 : lHanggi et all 120051 ). Among the necessary pre- 
requisites to run such heat machinery are thermal noise, 
nonlinearity, unbiased nonequilibrium driving of deter- 
ministic or stochastic nature and a symmetry break- 
ing mechanism. This then carries the setup away from 
thermal equilibrium, thereby circumventing the second 
law of thermodynamics, which otherwise would impose 
a van ishing directed transport ( Hanggi and Marchesonil 
12009(1 . 

Dwelling on similar ideas used in Brownian motors for 
directing particle flow, an efficient pumping or shuttling 
of energy across spatially extended nano-structures can 
be realized via modulating either one or more thermal 
bath temperatures, or applying external time-dependent 
fields, such as mechanical/electric/magnetic forces. This 
gives rise to intriguing phononic phenomena such as a 
priori directed shuttling of heat against an external ther- 
mal bias or the pumping of heat that is assisted by a 
geometrical-topological component . 



A. Classical heat shuttling 

In the following we shall elucidate the objective for 
shuttling heat against an externally applied thermal bias. 
A salient requirement for the modus operandi of heat mo- 
tors is the presence of a spatial or dynamic symmetry 
breaking. 

A possible scenario consists in coupling an asym- 
metric nonlinear structure to two baths; i.e., a left(L) 
and right (R) heat bath which can be modeled by 
classical Langevin dynamics. Applying a periodically 
time-varying temperature in one or both heat baths, 
Tl(r)(£) = Tl(r)(^ + 27r/a;), possessing the same aver- 
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FIG. 12 (color online). Action of a heat Brownian mo- 
tor in two coupled asymmetric FK-FK lattices. The heat 
baths are subjected to periodic modulations in the form 
T h (t) = T [l + A + A sgn(sinarf)] and T R (t) = T (l - A). 
The dimensionless reference temperature is set as To = 0.09 
(see the Appendix IA.2I for the expressions of corresponding 
dimensionless units). Note that in distinct contrast to Fig. [3] 
(d) the ratchet effect now yields with a modulation strength 
A 7^ a nonvanishing heat flow at zero bias A = 0. Apply- 
ing a substantial rocking strength, A = 0.5, the current bias 
characteristics can be manipulated as to inhibit a negative dif- 
ferential thermal resista nce (NDTR) re gime at negative bias 
values A. Adapted from lLi et al\ (|2008h . 



age temperature T L ( R )(i) — To, then brings the system 
out-of-equilibrium. Noteworthy is that this so driven sys- 
tem is unbiased; i.e., it exhibits a vanishing average ther- 
mal bias AT(t) = T L (t) - T R (t) = 0. The asymptotic 
heat flux J(t) assumes the periodicity of the external 
driving period 2ir/uj and the time-averaged heat flux J 
follows as the cycle-average over a full temporal period: 

J = /jf"" J(t)dt. Consequently, a nonzero average 
heat flux J ^ emerges which then provides the se ed to 
even shuttle heat against a net thermal bias AT(t) ^ 0. 

A first possibility to introduce the necessary symme- 
try breaking is to use an asymmetric material such as 
two coupled FK-FK lattices where two segments possess 
differ ent thermal properties ( Li et all 120081 : iRen and Lit 
I2010D . see in Fig. [T2J The directed heat transport 
can be extracted out of unbiased temperature fluctua- 
tions by harves ting the static thermal rectification effect 
(|Li et aLl . l200i . 

Yet another possibility is to break the symmetry dy- 
namically by exploiting the nonlinear response induced 
by the harmonic mixing mechanism, stemming from a 
time- varying two-mode modulation of the bath tempera- 
ture^), i.e., T L(R )(i) = T [l ± A\ cos(cj£) ± A 2 cos(2wi + 
tp)]. The second harmonic driving term A 2 cos (2a;f + tp) 
cause s the intended nonlinear frequency mixing (jLi et ail , 
120091) . 

In the adiabatic limit; i.e., if to — > 0, a nonzero heat 
flux J can be obtained due to the nonlinearity of 
the heat conductance response. The preferred direction 
of such heat flow is determined by the intrinsic thermal 
diode properties. In contrast, in the fast driving limit 
uj — > oo, the left- and right-end of the system will es- 
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sentially experience a time-averaged constant tempera- 
ture. This then mimics thermal equilibrium, yielding 
J — > 0. Remarkably, by modulating the driving fre- 
quency uj through the characteristic thermal response 
frequency of the system, the intriguing phenomenon of 
a hea t current reversal can be observed ( Li et all I2008L 
2009). For this two segment system, an optimal heat cur- 
rent can be obtained around this characteristic frequency 
even when the two isothermal baths are oscillating in 
synchrony with T L (t) = T R (i), and the cur rent reversa l 
can b e realized by tuning the system size ( Ren and Lit 
l2010l) . In the harmonic mixing mechanism the directed 
heat current is found to be proportional to the t hird-order 



mom ent (AT(t)/2T ) 3 , i.e., J oc A(A 2 cos ip (|Li et all 



I2009D . This enables a more efficient way of manipulating 
heat: the direction of heat current can be reversed by 
merely adjusting the relative phase shift tp of the second 
harmonic driving. 

Apart from using the FK-lattice as a source of nonlin- 
earity, other schemes of heat motors can be based on a 
Fermi-Pasta-Ulam (FPU) la ttice structure, a Lennard- 
Jones interaction pote ntial dLi et all 20091). or also a 
Morse lattice structure (|Gao and Zhengf 20 111) . 

Besides a manipulation of bath temperatures, the shut- 
tling of heat can be realized by use of a combination of 
time-dependent mechanical control in otherwise symmet- 
rical structures. Depending on specific nonlinear system 
setups, an emerging directed heat current can be con- 
trolled by adjus ting the relative pha se among the act- 
ing drive forces (IMarathe et "aZl. 120071) or the driving fre- 
quency (lAi et all 120101) . Multiple resonance structures 



for the heat current vs. t he driving frequency of external 
forces can occur as well (|Zhang et aU l2011bl) . It can be 
further demonstrated that for strict harmonic systems, 
periodic-force driving fails to shuttle heat. Moreover, it 
has been shown that even for anharmonic lattice seg- 
ments composed of an additional energy depot, it is not 
possible to pump heat from a cold reservo ir to a hot reser- 
voir by merely applying e xternal forces ()Marathe et all 
120071: IZhang et aLl . l2011bh . 

It is elucidative to compare these setups 
with the Feynman ratchet-and-pawl device 
of a heat pump (IVan den Broeck and Kawai , 



2006 



2008; 



.van den Brock and Van den Broeckl. 
Hang gi and Marchesoni . 20091 : Komatsu and Nakagawa , 
2006)~ The latter is a consequence of Onsager's recipro 



cal relation in the linear response regime: if a thermal 
bias generates a mechanical output, then an applied 
force will direct a heat flow as a conjugate behavior. 
Therefore, conjugated processes can be utilized for heat 
control as well. Other well-known thermally conjugated 
processes are the Seeback, Thomson and Peltier effects 
in thermoelectrical devices, where the thermal bias 
induces electrical currents, or vice versa. Recently, such 
a nanoscale magnetic heat engine and pump has been 
investigated for a magneto-mechanical system, which 
either operate as an engine via the application of a 
thermal bias to convert heat into useful work, or act as 
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FIG. 13 (color online). Quantum heat shuttling, (a) Setup 
of a molecular junction acting as a heat shuttle. The short 
molecular wire is composed of a single electronic level Ei only 
that can be gated and a single phonon mode at the fixed vi- 



brational frequency Wi = 1.4 x 10 



typical for a carbon- 



carbon bond. The lead temperatures r L ( R )(t) are subjected 
to time-periodic modulations, (b) Quantum thermal rachet 
effect for this molecular junction. Total directed heat current 
Jq as the function of static thermal bias ATbi as for different 
driving amplitude strengths A. The heat baths are modulated 
as T L {t) = T + \ AT bias + A cos Sit and Tn(t) =T - \ AT bias . 
The reference temperature is set as To = 300 K and the 
electronic wire level is set as Ei — n — 0.138 eV where p 
is the chemical potential. Jq is independent of the driving 
frequency Q as a cons equence of ballisti c heat transfer in the 
adiabatic limit. After IZhan et al\ (12009). 



a cooler via applying magnetic fields or mechanical force 
fields to pump heat (|Bauer et al .1. 120101) . 



B. Quantum heat shuttling 

The efficient shuttling of heat via a time-dependent 
modulation of bath temperatures can be extended into 
the quantum regime when tunneling and other quan- 
tum fluctuation effects come into play. In clear con- 
trast to the realm of electron shutt l ing dGalperin et al. . 
20071: lHanggi and Marchesonil 120091 lHanggi et all 12002 : 



Joachim and Ratnerl 120051 : iRemacle et all 2005), how 
ever, this aspect of shuttling quantum heat is presently 
still at an initial stage, although expected to undergo 
increasing future activity. 



1. Molecular wire setup 

In the following we consider one specific such case 
in some detail. Let us consider a setup with a typical 
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molecular wire for which the heat transport is generally 
governed by both, electrons and phonons. A schematic 
setup base d on a stylized mo lecular wire is sketched with 
Fig.fTgta) (|Zhan et aLl . l2009l ). The single electronic level 
Ei can be conveniently modulated by a gate voltage and 
wi denotes the vibrational frequency for a single phonon 
mode. The lead temperatures Tl(r)(^) undergo an adia- 
batically slow periodic modulation for both electron and 
phonon reservoirs. The latter is experiment ally acces- 
sible by use of a heating/cooling circulator ( Lee et all , 
2005). ' In an adiabatic driving regime, the asymptotic 

ballistic electron and phonon heat currents Jq < - ph ' ) (t) can 
be ca lculated via the celebrated Landauer- l ike ex pres- 
sions (jDubi and Pi Ventral . I2O1H : ISegal et~aD . l2003h : i.e., 



de 
2irh 



( £ - M )T el ( £ ) [/( £ ,T L (t))-/( £ ,T R (i))] 
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2tt 

where T cl (e) and T ph (u;), respectively, denote the tem- 
perature independent transmission probability for elec- 
trons with energy e and phonons with angular fre- 
quency to. Here, the functions /(e,T;(i)) = [exp((e — 
MOAeTlWRir 1 andn(w,Ti(«)) = [exp(fuv / k B Ti(t)) - 
1] _1 where / = L,R, denote the Fermi-Pirac distribution 
and the Bose-Einstein distribution, respectively. These 
functions both inherit a time-dependence which derives 
from the applied adiabatic periodic temperature modu- 
lation. 

Then, a finite total heat current Jq emerges, reading 



(8) 



This heat current Jq results as a consequence of the 
nonlinear dependence of quantum statistics on tempera- 
ture, see in Fig. [l~3l b): - note that in presence of modula- 
tion Jq is nonzero even for vanishing thermal bias ATbias ■ 
This ballistic heat transport is a pure quantum effect 
which will not occur in the classical diffusive limit, being 
approached at ultrahigh temperatures. The efficient ma- 
nipulation of heat shuttling can be realized by applying 
the above-mentioned harmonic mixing mechanism. Since 
the heat transport is carried also by electrons, adjusting 
the gate voltage gives rise to a intriguing control of heat 
current with the result that the direction of heat current 
experiences multiple reversals. 

The quantum heat shuttling of a dielectric molecu- 
lar wire can also be achieved upon periodically mod- 
ulating the molecular levels while this molecular wire 
is connected to two heat baths th at are character- 
ized by distinct spectral properties ( Segal and Nitzanl . 
2006). Interestingly, the pumping of quantum heat 
can be operated arbitrarily close to the Carnot effi- 
ciency by a tailored stochasti c modulation of the molec- 
ular levels (|Segall 120081 l2009h . Time-dependent phonon 
transport in the non-adiabatic regime and strong driv- 
ing perturbations has numerically been investigated by 
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FIG. 14 (color online). Quantum heat pumping via a ge- 
ometrical phase, (a) A schematic representation of a sin- 
gle molecular junction, (b) Quantum heat transfer across 

(g)the molecular junction is generated via an adiabatic two- 
parameter variation of the left bath temperature Tt(t) and 
the right bath temperature Tn(t), which maps onto a closed 

(7)(blue) circle. The arrow indicates the direction of the mod- 
ulation protocol, (c) Geometrical phase induced heat cur- 
rent , J gcom 5 versus the angular modulation frequency f2. The 
straight line corresponds to the analytic result while the open 
circles denote the simulation result s. For further details we 
refer the reader to iRen et al\ (|201fj ). 



use of t he nonequilibrium Gr e en's f u nction (NEGF) ap- 
proach (|Pubi and Pi Ventral 1201 it IWang et all 120081) 
for the case of coupled harmonic oscillator chains: 
There, the coupling to the reservoirs is held at differ- 
ent temperatures with the latter switched on suddenly 
( Cuansing and Wand . l2010h . 



2. Pumping heat via geometrical phase 

As discussed above, a one-parameter modulation, e.g. 
via the left-sided contact temperature Ti(t) is typically 
sufficient for quantum mechanical shuttling and rectifi- 
cation of heat, as sketched with Fig. H~4T a). 

It should be noted, however, that a cyclic modulation 
involving at least two control parameters generally in- 
duces additional current contributions beyond its mere 
dynamic component. This is so because with a varia- 
tion evolving in a (parameter) space of dimension d > 2 
one typically generates geometrical properties (i.e. a 
nonvanishing, gauge invariant curvature) in the higher- 
dimensional parameter space of the governing dynami- 
cal laws for the observable, in this case the heat flow. 
This geometrical properties in turn affect the resulting 
flow within an adiabatic or even non-adiabatic parame- 
ter variation. With a cyclic adiabatic variation of multi- 
ple parameters this contribution to the emerging flow is 
thus of topological origin. It has been popularized in the 
literature under the label of a geom etrical, finite c urva- 
ture or (Berry)-phase phenomenon ( Sinitsvnl l2009t ) . As 
a consequence, one needs to consider the total heat flow 
to be composed of two contributions, reading: 
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(9) 



Here, the geometrical contribution is proportional to 
the adiabatic small modulation frequency J7; implying 
that this correction is typically quite small in comparison 
with its non-vanishing dynamical contribution. There- 
fore, to observe this component it is advantageous to use 
a two-parameter variation such that the dynamic com- 
ponent vanishes identically. 

A sole geometrical contribution can be implemented, 
for example, by modulating the two bath temperatures 
in a way that the trajectory in the plane spanned by the 
two temperatures describes a circle, c.f. Fig. H3Tb). - in 
which case there results a sole Berry-phase-induced heat 
current , J gC om, while its dyna mical component , J r \ vn , is 
identically vanishing, Jd yn = ( Ren et al\ . I2010T ). Also, 
different from the case of an irreversible, dynamical heat 
flux, this geometrical contribution can be reversed upon 
simply reversing the protocol evolution. The latter op- 
eration thus provides a novel and convenient meth od for 
controlling energy transport. In iRen et all (|2010l ) the 
authors have demonstrated such nonvanishing quantum 
heat pumping as the result of a nonvanishing geometrical 
phase. This situation is depicted with Fig. [Iffi c). 

A similar geometrical phase effect can also be present 
in classical setups; e.g., for coupled harmonic oscillators 
in contact with Langevin heat baths. In this case it has 
been demonstrated that with a modulation of the two 
bath temperatures in time the geometrical phase phe- 
nomenon emerges only for the higher order moments of 
the heat flow; that is to say only beyond the average heat 
flux. Only when nonlinearity or temperature-dependent 
parameters in an interacting system are present can the 
geometrical phase manifest itself in producing a nonvan- 
ishing heat current. 

Moreover, the finite Berry-phase heat pump mecha- 
nism in both, quantum and classical systems has been 
demonstrated to cause a breakdown of the so termed 
"heat-flux fluctuation theorem", the latter being valid 
for a time-indepe ndent heat flux transfer. This fluc- 
tuati on theorem ( Campisi et all . 1201 it ISaito and Dharl . 
I2007D can be restored only under special co nditions in 
the pr esence of a vanishing Berry curvature IRen et all 
(j2010D . 



C. Topological phonon Hall effect 

It is known that a geometrical Berry phase yields pro- 
found effects on ele ctronic tran s port p roperties in various 
Hall effect setups ( Xiao et "all . 12010). Due to the very 
different nature of electrons and phonons, the phonon 
Hall effect (PHE) has b een discovered only recently in a 
paramagnetic dielectric ( Strohm et aZ.L [2005h . and subse- 
quently confirmed by yet a differ ent experimental setup 
( Invushkin and Taldenkovl 12007). In particular one ob- 
serves a transverse heat current in the direction perpen- 




FIG. 15 (color online). A schematic illustration of the 
phonon Hall effect with heat flowing along the direction indi- 
cated by the (red) arrow. 



dicular to the applied magnetic field and to the longitu- 
dinal temperature gradient, see in Fig. 1151 The discovery 
of this novel PHE renders the magnetic field to be an- 
other flexible degree of freedom for phonon manipulation 
towards the objective of energy and information control 
in phononics. 

Since the n, several theoretical e x planat io ns have been 
prop o sed dKaean and Maksimovl 2008 ;_ Sheng et~al\ . 



120061 : IWang and Zhang) . 120091 IZhang et all l2010d ) to 
understand the PHE by considering the spin-phonon 
coupling. The spin-phonon coupling has two possible 
origins: (i) It either derives from the magnetic vec- 
tor potential in ionic crystal lattices, where the vi- 
bration of atoms with effective charges wi ll experienc e 
the Lorentz force under magnetic fields ( Hold . Il972h . 
or (ii) it results from a Raman (spin-orbit) interaction 

iar 



(lloselevich and Capellmannl . ll99a[Kornig|.ll939l:IOrbachl 
119611 : IVan Vlecki . ll940D . It has been shown that by intro- 
ducing spin-phonon couplings, a ballistic system without 
nonlinear interact ion even exhibits the possibility of ther- 
mal rectification ( Zhang et aZ.L l2010bh . 



Similar to the various Hall effects occurring for elec- 
trons, a to pological exp l anation of the PHE has been 
provided in IZhang efaU (|2010aL l2011ah . The heat flow 
in the PHE is ascribed to two separate contributions: 
the normal flow responsible for the longitudinal phonon 
transport, and the anomalous flow manifesting itself as 
the Hall effect of the transverse phonon transport. A 
general expression for the transverse phonon Hall con- 
ductivity is obtained in terms of the Berry curvature 
of phononic band structures. The associated topologi- 
cal Chern number (a quantized integer) for each phonon 
band is defined via integrating the Berry curvature over 
the first Brillouin zone. For the two dimensional hon- 
eycomb and kagome-lattice, the authors observed phase 
transitions in the PHE, which correspond to the sudden 
change of the underlying band topology. The physical 
mechanism is r ooted in the touching and s plitting of the 
phonon bands (jZhang et aZ.I . l2010alT20TTah . 

Ther efore, much alike for ele ctrons in topological insu- 
lators ( Hasan and Kand . |2010|) , the design of a family of 
novel phononic devices - topological thermal insulators 
- is promising, with the bulk being an ordinary thermal 
insulator while the edge/surface constitutes an extraor- 
dinary thermal conductor. 
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V. SUMMARY, SUNDRIES AND OUTLOOK 

With this Colloquium we took the reader on a tour pre- 
senting the state of the art of the topic termed 'Phonon- 
ics', an emerging research direction which is expected 
to chime in the future with conventional 'Electronics'. 
Particularly, we surveyed and explained various physical 
mechanisms that are exploited in devising the elementary 
phononic toolkit: namely a thermal diode/rectifier, ther- 
mal transistor, thermal logic gate and thermal memory. 
These building blocks for doing phononic electronics are 
rooted in the application of suitable static or dynamic 
control schemes for shuttling heat. We further reviewed 
recent attempts to realize such phononic devices which 
all are based on nanostructures and discussed first exper- 
imental realizations. 



A. Challenges 

In spite of the rapid developments of phononics in both 
science and technology, we should stress that this objec- 
tive is still at its outset. To put these phononic devices to 
work, there remain lots of ambitious theoretical and se- 
vere experimental challenges to be overcome. For exam- 
ple, much work is still required for the physical realization 
of the phononic toolbox and for entering the next stage 
of assembling operating networks that are both scalable 
and stable under ambient conditions. 

Theory- Interface thermal resistance. As we have 
pointed out the underlying mechanism for thermal rec- 
tifier/diode is based on an asymmetric interface thermal 
resistance. However, thus far a truly comprehensive the- 
ory for this effect is lacking. 

The current approaches for thermal transport across 
an i nterface , such as the acoustic mismatch (AMM) the- 
ory (iLittld. Il959l) and the diff usive mismatch (DMM) 
theory ( Swartz and Pohll Il989h . are based on the as- 
sumption that phonon transport proceeds via a com- 
bination of either ballistic or diffusive transport on ei- 
ther side of the interface. Both schemes offer lim- 
ited accu racy for nano s cale i nterfacial resistance pre- 
dictions ( Stevens et all 120051) . due to the neglect of 
the atomic details of actual interfaces. Specifically, 
the acoustic mismatch model assumes that phonons 
are transported across the interface without being scat- 
tered; i.e. they are ballistic while the diffuse mismatch 
model assumes the opposite, namely that the phonons 
are scattered diffusively. Thus, the effects of scatter- 
ing on the interfacial thermal resistance act as upper 
and lower limits for real situation. In fact, both nu- 
merical and experime ntal studies in nanostructures rang- 
ing f rom nanotube s (IChang et ai , [2008; Zh ang and Lit 
2005), nan owires (lYang et all |2010|), to polyethylene 
nanofibers (|Henrv and Chenl . 120081 120091: IShen et all 
|2010[ ) all show that phonons undergo anomalous diffusion 
- i.e. so termed superdiffusion -, being faster than normal 
diffusion but slower than ballistic transport. Therefore, 



it is necessary to establish an improved theory describ- 
ing thermal transport across the interface by taking into 
account the anomalous thermal transport characteristics 
of nanostructures. 

Yet another fact of thermal rectification and neg- 
ative differential thermal resistance is the role of 
(strong) nonlinearity. How to set up a transport the- 
ory by incorporating nonlinearity in both, the quan- 
tum regime and the classical regime presents a chal- 
lenge. The non-equi l ibrium Green's (NEGF) function 
method ( Wang et all H008f) certainly serves as an ele- 
gant mathematical framework. However, when sophis- 
ticated phonon-phonon interactions (which derive from 
anharmonicity) , and alike, become increasingly impor- 
tant (this being so in realistic nanostructures), the NEGF 
presents a cumbersome task. In the classical regime, an 
effective phonon theor y has proven useful in character - 
izing heat conduction ( Li and Li . l2007t iLi et all 200Gb). 
Several cla s sical studies have been advanced (|He et all , 
I20101 . 120081 l2009h . which may as well be extended into 
the quantum regime. 

Experiment.-Eixjperimental realization of practical 
phononic device depends on how to measure its thermal 
conductivity accurately. To this end, one needs to get 
rid of the contact thermal resistance. Thus far no well 
defined scheme capable to eradicate such contact thermal 
resistance has been put forward. 

It also should be noted that presently the prime 
elementary phononic building block, i.e. the ther- 
mal r ectifier, has been realized on the micrometer 
scale (IChang et all 120061) and on the millimeter scale 
(|Kobavashi et all 120091: ISawaki et ai\ , 1201 ll ) only. The 
challenge for experimentalists is to make the sample 
smaller, for example towards a dozens of nanometers or 
even a few nanometers only. This then would boost the 
rectification to much larger values and simultaneously 
would allow detecting the asymmetric interface thermal 
resistance. The primary challenge ahead, however, is to 
validate the negative differential thermal resistance, i.e. 
the key element necessary in realizing the thermal tran- 
sistor. 

For electronic-like function a networking towards a 
phononic-like logic operation is indispensable. In the 
event that this task could be achieved successfully could 
then facilitate the desirable operation of an all-phononic 
computer. 



B. Future prospects 

As emphasized above, the area of phononics lingers 
still in its infancy, yet it is at the verge of blossoming up. 

From phononics to acoustics, and vice versa. - As 
we elucidated throughout this Colloquium the physical 
principle of phononic devices rests in the manipulation 
of phonon bands/spectra. This idea can be generalized 
to control any elastic/mechanical energy. For example, 
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inspired by the t hermal diode , an acoustic diode has 
been proposed by iLiang et all (2009) and subsequently 
realized experimentally by us e of a nonlinear a coustic 
medium and phononic cr ystal (iLiang et aLl.l2010|). and a 



sonic crystal geometry ( Li et al . 2011 ). iBoechler et all 



( 201lh have demonstrated experimentally elastic energy 
switching and rectification. There is little doubt that 
in parallel to phononics, namely an acoustic transistor, 
logic gate and possibly even a computer device may be 
realized in the foreseeable future. The ideas and con- 
cepts in controlling acoustic waves can as well be ex- 
ploited for phononics. For example, phononic crystals 
have been demonstrated t o be usefu l in m anipulating 
acoustic wave propagation ( Liu et all . l2000[ ) . This con- 
cept has been extended recently to control heat flow on 
the nanoscale (iHopkins efall. \201l\ lLandrv et all , 120081: 



iMcGaughev et all . 120061: lYu et aZ.I . l2010h . It would not 
come as a too big surprise when someday concepts such as 
a "heat-cloak" , "super heat lens" , etc., come alive via the 
implementation of phononic concepts by use of acoustic 
meta-materials (iFang et al. L 2006 ; lGuenneau et all . l2007t 
lYang et all . 12004 IZhang and Liul . 120041 ) . In facT active 
research is presently pursued in controlling and manipu- 
lating mechanical/elastic energy which ranges from long 
wavelength elastic waves to very short wavelength ther- 
mal waves, a glimpse of this activity can be adapted from 
the Proceedings of the wor l d first conference on 'Phonon- 
ics' (|E1-Kadv and Husseinl . l201ll) . 

PhoXonics: phononics plus photonics, and beyond. 
- Another promising prospect is based on the com- 
bination of phononics with photonics to control and 
manage photon and phono n energy concomitantly 
(|Maldovan and Thomasl . 120061) . This synergy might po- 
tentially enable one to use the solar energy more wisely. 
For instance, a single-molecule phonon field-effect tran- 
sistor has been designed wherein phonon conductance is 
contr olled by a back- gate electric field ( Menezes et all . 
|2010| ). Moreover, as temperature is the commonly ap- 
plied control parameter for chemical and biological reac- 
tions, phononic devices may also find application for the 
local control of temperatures, for example in regulating 
molecular self-assembling processes with the "Lab-on-a- 
Chip" technology. 

Phononics and electronics. - Electrons carry as well 
heat. Therefore one can control the heat flow carried 
by electrons with the help of electric/magnetic fields. 
For example, by asymmetrically coupling a quantum 
dot to its two leads o ne can construct a heat rectifier 
( Scheibner et all . 120081) : - upon applying a gate voltage 
it is possible to operate a heat transistor ( Saira et all . 
120071 ). These kind of devices combined with phononic 
circuits may then carry the potential to manipulate dis- 
sipation of heat and cooling in nanoscale and molecular 
devices. Overall, there is also the potential that hybrid 
structures composed of electronic and phononic elements 
may lead to beneficial applications. 

From heat conduction to radiation.- In this colloquium, 
our proposed phononic devices are based on the control 



of heat conduction, being assisted by lattice vibrations. 
A similar idea can be generalized to control heat radia- 
tion. Indeed, Fan and collaborators at Stanford proposed 
a photon-m e diated "thermal rectifier through vacuum" 
(lOtev et all |2010D which makes constructive use of the 
temperature dependence of underlying electr o- magnetic 
reson ances. In the same spirit, Fan's group ( Zhu et "all . 
20121 ) revealed negative differential thermal conductance 
through vacuum; this in turn allows for the blueprint of 
a transistor for heat radiation. 

Finally, let us end with a celebrated quote by Winston 
Churchill: "This is not the end. It is not even the be- 
ginning of the end. But it is, perhaps, the end of the 
beginning." 
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Appendix: Nonlinear Lattice Models 
1. Lattice models 

In this Appendix we introduce three archetype one- 
dimensional (ID) lattice models commonly used in the 
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investigation of heat transport. These are (i) the linear 
Harmonic lattice, (ii) the nonlinear Fermi-Pasta-Ulam j3 
(FPU-/?) lattice and (hi) the Frenkel-Kontorova (FK) lat- 
tice. The simplest harmonic lattice serves as the basic 
model from which the FPU-/? lattice and FK lattice are 
derived by complementing the dynamics with a nonlin- 
ear inter-atom interaction in the FPU-case and an on-site 
potential in the FK-case, respectively. 

For a ID harmonic lattice with N atoms the normal 
modes of the lattice vibrations are known as phonons. 
The corresponding harmonic lattice Hamiltonian explic- 
itly reads 



reading: 



/v 



* = E 



Pt 
2m 



fc 



(Al) 



wherein the dynamical variables pi and x% i — 1,...,N, 
denote the momentum and position degrees of freedom 
for the i-th atom and xq = x\ — a. The parameters to, 
ko, a denote the mass of the atom, the spring constant 
and the lattice constant, respectively. The position vari- 
able Xi can be replaced by the displacement from equi- 
librium position as Sxi = xi — ia which we denote by the 
same symbol Xj = Sxi henceforth. The Hamiltonian of 
Eq. (|A1[) with Sxq = 5x\ thus simplifies, reading: 



JV 
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(A2) 



This form explicitly depicts the translational invariance 
of the free chain which implies momentum conservation. 

Applying next periodic boundary conditions xq = xn, 
the harmonic lattice of Eq. (|A2|) can be decomposed into 
the sum of noninteracting normal modes (phonons) with 
the dispersion relation reading 



u(q) = 2\ — |sin(g/2)| , < q < 2ir (A3) 

TO, 



where the continuous spectrum is due to the adoption of 
thermodynamical limit N —¥ oo. The harmonic lattice 
possesses an acoustic phonon branch with ui(q) —> as 
q -> 0. 

Although the 3D harmonic lattice model yields a sat- 
isfactory explanation for the temperature dependence of 
experimentally measured specific heat, it turns out that 
this very model of noninteracting phonon modes fails to 
descri be a Fourier Law for heat transport. In pioneering 
work ( Rieder et all Il967l ) the authors proved that for 
this ID harmonic lattice the heat transport is ballistic; 
this is due to the absence of phonon-phonon interactions. 
In order to take the phonon-phonon interactions into ac- 
count, the harmonic chain model must be complemented 
with nonlinearity. 

If a quartic inter-atom potential is added, one arrives 
at the FPU-/? lattice with the corresponding Hamiltonian 
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where the parameter j3o is the nonlinear coupling 
strength. Historically, the FPU-/3 latti ce has been 
put f orward by Fermi, Pasta and Ulam ( Fermi et all 
1955) to study the issue of ergodicity of a nonlin- 
ear system dynamics. For some excellent reviews of 
the o riginal FPU problem we refer the readers to 
Refs. (d crman and IzrailevL 120051: iFordl Il992ri ). Surpris- 
ingly, the FPU-/? lattice still fails to obey Fourier's law 
and the heat conductivity k diverges with th e system size 
propo rtional to k k N a with < a < 1 (|Lepri et all 
Il997f ) . This divergent behavior has recently been experi- 
mentall y verified for a syst em setup using quasi- ID nan- 
otubes |chang et aEl . l2008| ). The phonon modes of FPU- 
/? lattice are also acoustic-like after re-normalization of 
the nonlinear part, due to the conservation of total mo- 
mentum. 

If a periodic on-site substrate potential is added to 
the harmonic chain, one arrives at the FK lattice. Its 
Hamiltonian is given by 
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where parameter Vq/Att 2 denotes the nonlinear on-site 
coupling strength. Here we only consider the commen- 
surate case where the on-site potential assumes the same 
spatial periodicity as the harmonic lattice. Notably this 
model with an on-site potential now breaks momentum 
conservation. Among the various phenomenological mod- 
els that mimic solid state systems the FK model has 
been shown to provide a suitable theoretical descrip- 
tion for possible nonlinear phenomena such as the oc- 
currenc e of commensura t e-inco mmensurate phase tran- 
sitions ijFlori^me^Maz^!. 19961). kink- like structures and 



alike (jBraun and Kivshaii Il998l . [20041 ). It has attracted 
interest since it was first proposed by Frenk el and Kon- 
torova (|Frenkel and Kontorovai Il938l H939) in order to 



study various surface phenomena. Recently it has been 
established that the FK lattice indeed does exhibit nor- 
mal heat condu ction and thus obeys the Fourier law 
( Hu et all Il998f) . This normal behavior is attributed to 
the optical phonon mode where the phonon mode opens 
a gap as the momentum conservation is broken with the 
on-site potential. 



2. Local temperature and heat flow 

Dimensionless units constitute practical tools for the 
theoretical analysis and numerical simulations. Here we 
provide a brief introduction to the dimensionless units 
used in this report for the various lattice model setups. 
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Let us start with the simplest lattice model of ID Har- 
monic lattice of Eq. (|A2[) . For the harmonic lattice con- 
tacting a heat bath specified by a temperature T, there 
are four independent parameters m, a, ko and fee where 
denotes the Boltzmann constant. The dimensions of 
all the physical quantities that typically enter the issue 
of heat transport can be expressed by the proper com- 
bination of these four independent parameters because 
there are only four fundamental physical units involved: 
length, time, mass and temperature. 

As a result, one can introduce the dimensionless vari- 
ables by measuring lengths in units of [a], momenta in 
units of [a(mfco) 1 / 2 ], temperature in units of [fcoaVfce], 
frequencies in units of [(fco/m) 1 / 2 ], energies in units of 
[k a 2 ] and heat currents in units of [a 2 k^ 2 /(27rm 1 / 2 )]. 
In particular, the Hamiltonian of Eq. (IA2[) can be 
transformed into a dimensionless form if we implement 
the following substitutions: 

H -> H[k a 2 ],pi -> pi[a(mk ) 1/2 },Xi -> Xi[a], (A6) 

where the so transformed dynamical variables yield the 
dimensionless variables to obtain 



N r 



1 



(A7) 



Typical physical values for atom chains are as follows: 
a ~ 10" 10 m, u} ~ 10 13 sec" 1 , m ~ 10~ 26 - 10~ 27 
kg, k B = 1.38 x 10" 23 J K _1 , we have [k a 2 /k B ] ~ 
(10 2 — 10 3 )K. This in turn implies that room temper- 
ature corresponds to a dimension less temperature T of 
the order 0.1-1 (|Hu et all Il998t l . 

To obtain the dimensionless FPU-/3 lattice from 
Eq. (|A4[) . one cannot scale the five parameters k B = o, = 
m = ko = 0o = 1 because one of them is redundant. 
Applying the substitutions of Eq. (|A6|) . we obtain the 
dimensionless form for the FPU-/? Hamiltonian: 
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with the dimensionless parameter (3 = /?oa 2 /fco- It is ev- 
ident that the dimensionless nonlinear coupling strength 
/3 is generally not equal to unity. However, it can be 
shown that upon adjusting /3 becomes equivalent to vary 
the system energy or its temperature. 

The dimensionless FK Hamiltonian can also be ob- 
tained by use of Eq. (IA6|) : 
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where the dimensionless on-site coupling strength V = 
Vo/koa 2 . 

Thus far, we dealt with homogeneous lattice Hamil- 
tonians. For thermal devices with more than one seg- 
ment, each segment may possess its own set of parameters 



such as a different spring constant or nonlinear coupling 
strength. In those cases, the reference parameter, for in- 
stance ko, used to define a transformation in Eq. (|A6j) 
may be chosen to correspond to a natural parameter of 
the corresponding segment. In particular, the dimension- 
less Hamiltonian for each individual segment of a coupled 
FK-FK lattice may be written as 
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where k is measured with the reference to a parameter 
ko which is introduced a priori. 

Next we discuss the results for expressing temperature 
and heat current in dimensionless units. In our classi- 
cal simulations we typically used Langevin thermostats 
with coupling the "contact" or end particles to the heat 
baths at their corresponding temperature. More pre- 
cisely, one adds to the corresponding Newtonian equation 
of motion a Langevin fluctuating term whic h satisfies the 
fluctu ation-dissipation relation, e.g. see in lHanggi et all 
(1990). Towards this goal we made use of the equipar- 
tition theorem of classical statistical mechanics to de- 
fine, for example, the local temperature Tj via its average 
atomic kinetic energy; i.e., 



fcT, = (£■)_> Ti = (pf) 



(All) 



where the arrow indicates the dimensionless substitu- 
tion Ti — > Ti[koa 2 /k B ] and pi — > pi[a(mfco) 1 ^ 2 ], and 
(• • •) denotes the (long) time-average or, equivalently, its 
ensemble-average in numerical simulations, thus implic- 
itly assuming ergodicity (in mean value) . 

Unlike temperature, the expression for the heat current 
is model dependent. To arrive at a compact expression 
we first rewrite the ID lattice Hamiltonian in the more 
general form: 



" = E 



(A12) 



where V(xi—±,Xi) denotes the inter-atom potential and 
U(xi ) is the on-site potential. As a result of the continu- 
ity equation for local energy, the local, momen tary heat 
current can be expressed as ()Lepri et a/1 l2003h : 



dV(xi-i,Xj) 
ox l 



(A13) 



Consequently, the expression of heat current depends 
only on the form of the inter-atom potential V(xj_i, Xj). 
One should notice that although the on-site potential 
U(xi) does not enter into the expression of heat current 
explicitly, it does, however, influence the heat current 
implicitly through the dynamical equations of motions. 

The heat currents itself are again obtained via the 
time average over an extended time span. For steady 
state setups with fixed bath temperatures the resulting 
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heat currents are time-independent and, as well, inde- 
pendent of the particular site index (i) within the partic- 
ular chain segment. Likewise, with periodically varying 
bath temperatures T(t), the resulting ensemble average 
is also time-periodic; an additional time-average over the 
temporal period of T(t) yields the cycle-averaged, time- 
independent heat flux. Alternatively, an explicit long 
time average again produces this very time-independent 
value for the heat current. 



3. Power spectra of FPU-/3 and FK lattices 

The power spectrum (or power spectral density) de- 
scribes the distribution of a system's kinetic energy 
falling within given frequency intervals. For a homoge- 
neous lattice composed of identical particles the velocity 
Vi(t) = v(t) of a particle located at site i becomes inde- 
pendent of position i; the power spectrum then can be 
conveniently calculated by the Fourier transform of the 
corresponding velocity degree of freedom to yield: 



lim — / ° v{t)e- lult dt 

to-J-oo t J 



(A14) 



In doing so, the power spectrum of the FPU-/3 model 
of Eq. (|A8p depends on the temperature. In the low 
temperature regime the FPU-/3 dynamics is close to a 
harmonic lattice, yielding < u < 2. In contrast, in the 
high temperature regime it is the anharmonic part that 
starts to dominate. In this latter regime an approximate 
theoretical estimate then yields < oj < (^(T/?) 1 / 4 , with 
C = 272^r(3/4)3 1 / 4 /r(l/4) » 2.23 where V denotes 
the Gamma function (|Li et "o?l l2005a|) . Therefore, upon 
increasing the temperature then causes a rightward shift 
of the power spectrum towards higher frequencies. The 
Parseval's theorem then dictates that the area below the 
curve is proportional to the average kinetic energy of 
the particle; i.e., J Q P{tu)du> ~ (Eki n )- 



For the FK model in Eq. (|A10[> the form of the power 
spectrum depends sensitively on temperature. In the low 
temperature limit, the atoms are confined in the valley 
of the on-site potential. Upon linearizing Eq. (|A10I) the 
phonon band can be extracted to read: 



V < uj < VVTIk. 



(A15) 



In contrast, in the high temperature limit the on-site 
potential can be neglected; thus the system dynamics be- 
comes effectively reduced to a harmonic chain dynamics, 
whose phonon band extends to: 



< u < 2\fk. 



(A16) 



The crossover temperature T c can be approximated as: 
T cr ps V/(2ir) 2 . Its value depends on the height of on- 
site potential. This in turn implies different values for 
the two segments of the thermal diode setup, see Fig. [121 
This difference is at the heart of the thermal rectifying 
mechanism. 




k=0.2 
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T=0.02 
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FIG. 16 (color online). Temperature dependent power spec- 
tra. The variation of the power spectrum at different tem- 
peratures vs. the angular frequency uj (both in dimensionless 
units) in an FK-lattice with 10000 sites for two different sets 
of parameters in (a) and (b) . The features of these nonlinear 
FK-power spectra provide the seed for the modus operandi in 
a thermal diode setup as discussed in Sec III.ATT1 The panel 
(a) is for a FK-coupling strength of V=5 and a strength for 
the spring constant of k=l; panel (b) is for a coupling strength 
V=l and a spring constant value set at k=0.2. 
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